Numerical study of scars in a chaotic billiard 
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Abstract 



a^ 

We study numerically the scaling properties of scars in stadium billiard. 
Using the semiclassical criterion, we have searched systematically the scars of 
the same type through a very wide range, from ground state to as high as the 
1 millionth state. We have analyzed the integrated probability density along 
the periodic orbit. The numerical results confirm that the average intensity 
O . of certain types of scars is independent of h rather than scales with Vh. Our 

I findings confirm the theoretical predictions of Robnik (1989). 
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Eigenstates of a bound quantum system whose classical counterpart is chaotic are of 
great interest in the fast developing field of quantum chaos. Among many others, scars are 
one of the most interesting and striking topics. Since its discovery [|1],0], much progress has 
been achieved. On the theoretical side, Bogomolny [|| developed semiclassical theory of 
scars in configuration space, and Berry performed a similar analysis in phase space using 
the Wigner function. According to this theory, the intensity of a scar goes as \fTi. Based 
on the semiclassical evaluation of the Green function of the Schrodinger equation in terms 
of the classical orbit, Robnik |^ has developed a theory, and predicted that if the scar is 
supported by many periodic orbits, the maximal intensity of the scar is independent of ^, 
although its geometry can be determined by Bogomolny's theory. Most recently, Klakow 
and Smilansky used a scattering quantization approach for this problem. Parallel to 
the theoretical developments, there have also been many numerical and experimental 
studies @. 

Unfortunately, due to the limit of the numerical techniques and the computer facilities, 
most of the numerical studies up to now are restricted to a very low energy range, which is 
too low to verify the theoretical predictions in the very far semiclassical limit, particularly 
for Robnik's theory. In this paper, by using our numerical code of the improved plane wave 
decomposition method (PWDM) (for more details about Heller's PWDM, please see [ITU 



while for the improved PWDM, we will discuss it in another paper [0), we are successful 
to go as high as the 1 millionth state, which is very deep in the semiclassical regime for the 
stadium billiard. With the help of the semiclassical criterion |P , we found many consecutive 
scars in several different energy ranges, which spans 2 orders of magnitude in the wave 
vector, therefore, we are able to study the properties of scars, such as the intensity and 
profiles in the very far semiclassical limit. To the best of our knowledge, this is the first one 
of such kind of study. 

To make the numerical data significant, we need enough ensembles of scars of the same 
type. Therefore, our first step is to collect scars of the same type in a wide range of energy. 
We begin from a very low state, e.g. the ground state. As long as we find the first scar, say, 
e.g., at wave vector fco, then we can use the semiclassical criterion to estimate the next scar. 
According to the semiclassical theory [H-H], the scar is most likely to occur if quantized, i.e., 



S = 2Txri{n^^, n = 0,1,2,.... (l) 

S is the action along the periodic orbit, a is the Maslov phase. It must be pointed out 
that the semiclassical theory cannot predict the individual state at which the scar will occur. 
Instead, if we have already found one scar, say, at fco, then the semiclassical theory tells us 
that the eigenstates at the wave vector of /cq ± A/c will be scarred most likely. AA; = 27r/i/£, 
L is the length of the periodic orbit. In our study we put ^ = 1, so the inverse of the wave 
vector k plays the role of ^, i.e., k goes to infinity indicates the semiclassical limit. It has 
been verified in our numerical study that this criterion is very helpful and very successful 
in searching for and collecting scars. As we shall see later in many cases this criterion is 
accurate within one mean level spacing, namely, the scar occurs at the eigenstate whose 
eigenenergy roughly equals the predicted energy. 

With the help of the semiclassical quantization criterion equation (1), we have found 
about 100 examples of the same type of scar at different energy ranges. One such example is 
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shown in Fig. 1. The eigenvalue of this eigenstate is k = 1328.153 849, which corresponds to 
the sequential number 250 034 for odd-odd parity, and to the index of about 1 001 408 when 
all parities are taken into account. To our surprise, in addition to this one, we have found 
quite a few examples of this type of scarred state in such a high level. This implies that 
the scars survive the semi classical limit. Does this finding contradict with the Shnirlman's 



theorem |Tl| , which states that as the energy goes to infinity, the probability density of most 
eigenstates of a chaotic billiard approaches a uniform distribution? To test this, we have 
investigated the statistics of the probability density distribution of the wave function, and 
found that it is an excellent Gaussian distribution, although there is a pronounced density 
around the periodic orbit. 

In order to understand the scar properties quantitatively, we have investigated the fol- 
lowing pronounced (excess) intensity in a thin tube along the periodic orbit (see Fig. 2), 
which is defined by 

, /z^2(x)rfx 



1. (2) 



where ip{x.) is the eigenfunction at x. (■?/'^(x)) is the average probability density inside 
billiard, which is 1/^ according to the semiclassical theory |Tl]-|13|. A is area of the billiard. 
The integral is taken over a thin tube around the periodic orbit, which is presented in Fig. 
2. 

According to Robnik's theory [Q, although the geometry of a scar is determined by 
a single short periodic orbit, the intensity profile is nevertheless determined by the sum 
of contributions from similar but longer periodic orbits, which "live" in the homoclinic 
neighborhood close to the stable and unstable manifolds of the primitive orbit. Taking 
into account all these orbits, the pronounced intensity of the scar defined by Eq.(0) can be 
described by the following formula, 

sin{nS,/h) ^ 
^^sinh{nXT/2 ' ^ ' 

where, 5*1 is the action along the primitive periodic orbit, A is the Lyapunov exponent of 
the primitive orbit with the period of r, the summation over n is due to the repetitions of 
the orbit, and p is the number of contribution orbits. Equation (|^) states that the maximal 
intensity of the scar, when supported by many periodic orbits, is independent of fi. This 
thoeretical prediction is different from that of Bogomolny. But, it does not contradict that 
of Bogomolny at all, instead, it is an extension of Bogomolny's theory to the scars caused 
by many periodic orbits. These two theories describe different types of scars. Indeed, we 
have also found the scar types whose intensity depends on a/^, which is exactly predicted by 
Bogomolny's theory. However, since other authors |p have already verified this theoretical 
prediction, we will not repeat this in this paper, we shall concentrate on the scars that 
cannot be determined by Bogomolny's theory, but can be described by Robnik's theory. 

In Fig. 3, we show six representative examples of this intensity versus the width of the 
tube {D) in units of de Broglie wavelength around the periodic orbit. These six examples 
are the same type of scar, namely, the diamond-shaped scar shown in Fig. 1. They go from 
the very low state k= 10.241 095 to the very deep semiclassical regime at /c = 1328.153 849. 
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The first thing one can see from these profile figures is that the scar intensity has a 
maximum at the width of about 1-2 de Brodglie wavelengths from the periodic orbit. This 
can be explained by Robnik's theory. The semiclassical waves associated with individual 
daughter orbits interfere constructively with each other only within a tube of width 1-2 de 
Brodglie wavelengths. The second important result from this figure is that the magnitude of 
the maximum does not change too much although the eigenenergy changes more than 100 
times. 

Moreover, after checking the eigenenergies of these six examples carefully, we have found 
that the semiclassical criterion works very well, although we go from one scar state to another 
by jumping even up to a few hundred scarred states. For instance, starting from the first 
eigenvector ko = 10.241 095, if we go through 65 scarred states, we have k = ko + 65Ak = 
101.563 684, which is very close to the exact one kexact =101.568 640. ( In this paper, we 
study only the eigenstates with odd-odd parity, so the length of the periodic orbit shown 
in Fig. 2 is £ = 2a/5, rather than 4^5 for the total biUiard, thus, Ak = 27r/£ = 1.404 
96). The deviation is less than one mean level spacing. This procedure applies also to many 
other scarred states and it can be verified readily for other states given in Fig. 3. 

The next very important question is that how does this maximal integrated intensity 
depend on the eigenenergy or the h7 Firstly, our numerical results show that around a 
certain k, it changes from scarred state to state. This is shown in Fig. 4, where we plot 
26 consecutive scarred states around k = 125. (Note that there are two cases in which 
two consecutive eigenstates are near degenerate, thus both of them are scarred.) Again, 
from this figure we can see clearly that the semiclassical criterion (|l|) works very excellently. 
The interval between two scarred states is almost constant and approximately equals 2n/C 
The maximal integrated intensity, however, fluctuates from state to state, which cannot be 
explained by the existing semiclassical approaches. This is still an open problem deserving 
of further theoretical and numerical investigations. 

The results shown in Fig. 4 imply that in order to make the study of the dependence 
of the maximal integrated intensity on energy significant we should take certain kinds of 
ensemble averaging. In our numerical study, we have performed such averaging around 
certain k over many scarred states (usually about 10 states). The averaged data are plotted 
in Fig. 5. The best least-squares fit gives rise to 

{Im) = 0.73/r, a = 0.06 ± 0.03, (4) 

here, () means the local average. The exponent a is very close to zero and is far from 1/2 
as predicted by Bogomolny's theory. This fact means that the maximal integrated intensity 
does not depend on the energy or h for the scar type shown in this paper. This discovery is 
very different from previous ones and cannot be explained by the semiclassical theory of 
Bogomolny ||^ and Berry 0, however, it confirms quantitatively the theoretical prediction 
of Robnik f^, which states that the maximal intensity of a scar is independent of h if the 
scar is supported by many orbits, as mentioned above. 

In this paper, we have studied intensively the scars in a stadium billiard, and have shown 
numerically that the semiclassical criterion (|I]) works very well from very low state to that 
in the very far semiclassical limit. Furthermore, we have analyzed the scaling property of a 
scar with h and found that for the scar type shown in this paper, the maximal integrated 
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density fluctuates from scarred state to state, but the local average intensity does not change 
with energy. This flnding conflrms Robnik's scar theory of multiple periodic orbits 

The author would like to thank Professor Dr. Marko Robnik for discussions. He is also 
very grateful to Professor Dr. Felix Izrailev for helpful discussions during the STATPHYS 
19 in Xiamen (1995) and during his visit in Como (1996). This work was supported in part 
by the Research Grant Council Grant RGC/96-97/10 and the Hong Kong Baptist University 
Faculty Research Grants FRG/95-96/II-09 and FRG/95-96/II-92. The work done in Slovenia 
was supported by the Ministry of Science and Technology of Republic of Slovenia. 



5 



REFERENCES 



[1] S.W. McDonald and A.N. Kaufman, Phys. Rev. Lett. 42, 1189 (1979); S.W. McDonald, 

Ph.D. Thesis, Berkeley, 1981. 
[2] E.J. Heller, Phys. Rev. Lett. 53, 1515 (1984). 
[3] E.B. Bogomolny, Physica D 31, 169 (1988). 
[4] M. Berry, Proc. R. Soc. London A 423, 219 (1989). 

[5] M. Robnik, Preprint of ITP, University of California at Santa Barbara, 1989. 
[6] D. Klakow and U. Smilansky, J. Phys. A: Math. Gen. 29, 3213 (1996). 
[7] R.L. Waterland, J.-M. Yuan, C.C. Martens, R.E. GiUilan, and W.P. Reinhardt, Phys. 
Rev. Lett. 61, 2733 (1988); B. Eckhard, G. Hose, and E. Polak, Phys. Rev. A 39, 3776 

(1989). 

[8] D. Wintgen and A. Honig, Phys. Rev. Lett. 63, 1467 (1989); O. Agam and S. Fishman, 
ibid 73, 806 (1994). 

[9] S. Sridhar and E. J. Heller, Phys. Rev. A 46, 1728 (1992); J. Stein and H.-J. Stockman, 
Phys. Rev. Lett. 68, 2867 (1992). 

[10] E.J. Heller, in Proceedings of the 1989 Les Houches Summer School on "Chaos and 
Quantum Physics", edited by M.J.Giannoni, A.Voros and J. Zinn- Justin (Elsevier Sci- 
ence Publisher B.V., Amsterdam, 1991), p602; B. Li and M. Robnik, J. Phys. A: Math. 
Gen. 27, 5509 (1994). 

[11] A.L. Shnirelman, Usp. Mat. Nauk. 29, 181 (1974). 

[12] M. Berry, J. Phys. A: Math. Gen. 12, 2083 (1977). 

[13] A Voros, in Stochastic Behavior in Classical and Quantum Hamiltonian Systems, edited 
by G. Casati and J. Ford; Lecture Notes in Physics 93, (Springer, Berhn 1979), p. 326. 

[14] B. Hu and B. Li, Statistical analysis of scars in stadium billiard, Preprint CNS-9703, 
Centre for Nonlinear Studies, Hong Kong Baptist University, 1997 



6 



FIGURES 



FIG. 1. The probability density plot for a scarred eigenstate of odd-odd parity. The wave 
number k = 1328.153 849, which corresponds to index 250 034 using the Weyl formula (odd-odd), 
thus it corresponds to approximately the 1 001 408th eigenstate for the total billiard. The scar is 
obviously supported by the diamond-shaped periodic orbit shown in Fig. 2. The stadium has the 
parameter of circle radius R = 1 and the straight line length 2. In this figure, the unit length is 
about 211 de Broglie wavelengths. 

FIG. 2. The integral region around the periodic orbit that is taken in Eq. (2). The width of 
the tube is D measured perpendicular to the periodic orbit. 

FIG. 3. The integrated scar intensity / vs the width of the integrating tube in unit of the de 
Broglie wavelength for the scar type shown in Fig. 1. We show six different eigenstates at different 
energy ranges. 

FIG. 4. The maximum of integrated scar intensity vs wave number k around k = 125. The 
type of scar is the same as shown in Fig. 1 (the diamond shape). Here we see clearly that the wave 
number interval between two consecutive scars is very close to 2tt/C (= 1.404 96), as predicted by 
the semiclassical quantization condition Eq.(|l|). 

FIG. 5. The locally averaged (over a small group of consecutive scarred states) maximum of 
integrated scar excess intensity vs wave number k. The bullet represents the numerical data, and 
the solid line is the best least-squares fit. 
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